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Abstract. The cyclic oxidation process can be described as an iterative scale growth and spallation 
sequence by a number of similar models. Model input variables include oxide scale type and 
growth parameters, spalling geometry, spall constant, and cycle duration. Outputs include net 
weight change, the amounts of retained and spalled oxide, the total oxygen and metal consumed, 
and the terminal rates of weight loss and metal consumption. All models and their variations 
produce a number of similar characteristic features. In general, spalling and material consumption 
increase to a steady state rate, at which point the retained scale approaches a constant and the rate of 
weight loss becomes linear. For one model, this regularity was demonstrated as dimensionless, 
universal expressions, obtained by normalizing the variables by critical performance factors. These 
insights were enabled through the use of the COSP for Windows cyclic oxidation spalling program. 

Introduction 

For isothermal exposures, degradation from oxidation is easily defined on the basis of scale growth 
models alone. However, in cyclic exposures, protective scales may be especially compromised 
because some portion of the scale may spall off on cooldown due to extremely high compressive 
thermal expansion misfit stresses. The primary effect is that, after portions of the scale are 
removed, re-growth of a thinner scale occurs at rates higher than those for continued growth of the 
undamaged thicker scale. This spalling amount and configuration is therefore just as important as 
growth rate in determining metal consumption. Rather than fit experimental weight change curves, 
the purpose of this paper is to highlight and review the typical trends predicted for oxidation and 
metal consumption as a function of oxide type, growth kinetics, spalling formalisms, and cycle 
duration. Universal response curves for one particular model are presented by using dimensionless, 
normalized parameters. 

Elements of Cyclic Oxidation Models 

Various models have been put forth to describe these processes and successfully fit experimental 
cyclic oxidation curves to model predictions [1-3]. Perhaps the most all-inclusive treatise is the 
COSP Cyclic Oxidation Spalling Program [2-4], This model applies to the majority of the topics 
discussed below, except where specifically noted. 

Scale Growth Models. Scale growth is normally characterized in isothermal tests, and parabolic 
rates are normally encountered for protective scales due to diffusion control of the growth process. 
Modifications due to scale grain growth coupled with grain boundary diffusion control or due to 
rapid transient oxidation occurring before protective healing layers are formed may lead to power 
law or logarithmic rates as listed below [2]: 


parabolic: (AW/A) 2 = (k p t) 

(1) 

power law: (AW/A) m = (kt) 

(2) 

logarithmic (AW/A) = In [(kt+c) 1/m ] 

(3) 
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Spalling Models. The geometry or “fractopography’ of the spalled portion must be specified in 
order to completely define the thickness of the retained scale and the corresponding growth rate 
upon reheating. At one extreme, all spallation is assumed to occur uniformly as an outer layer of 
the scale, Figure 1, Case A. This would apply to extremely adherent scales that remain attached to 
the substrate, then fracture within the scale after accumulating a sufficient amount of stored strain 
energy. While a single uniform layer is not observed in practice, it represents a simplifying 
assumption and approximates the response, on average, over the entire sample surface. 

At the other extreme is the use of discrete segments which do spall completely to the scale-metal 
interface, Fig. 1, Case C. And an intermediate situation, Case B, allows for segment spallation at 
some depth fraction of the scale thickness, R s . Finally a bimodal combination of A and B (or C), 
each having different depth fractions of spallation, R s , and each with a representative spalling 
probability, Pi and P 2 , is proposed. For these configurations, rather than picking the same segments 
to spall on every cycle, the occurrence is randomized over all segments by means of a Monte Carlo 
algorithm [2]. 



Figure 1. Various schemes of scale 
‘fractopography’ for COSP models. 

A) uniform surface layer; 

B) partial depth fraction, R s ; 

C) full interfacial segment; 

A + (B or C) bimodal, with 
probabilities of Pi and P 2 . 


Spalling For malisms . One of the earliest models employed a constant interfacial spalling area 
fraction, k s , wherein each successive cycle caused every segment to spall ‘k s ’ of its area, leaving an 
unspalled area fraction equal to (l-k s ) [1]. This model required the unusual morphology of having 
2 1 segments, where i refers to number of cycles, resulting in a relative area of (l-ks ) 1 of originally 
intact oxide. A more recent interfacial model, DICOSM, (Deterministic Interfacial Cyclic 
Oxidation Spalling Program), also specifies a constant area fraction, Fa, to spall at the interface, but 
only at the thickest portion of the scale [5]. But here the number of effective area segments is now 
fixed, as 1/Fa. A similar interfacial model has defined the spalling process by a statistical 
probability factor, p, which relates the overall tendency for the scale to spall on any given cycle 
[6,7]. Identical behavior is observed between the probabilistic, p, and the area fraction interfacial 
model, k s , but divergence occurs between these and the DICOSM, Fa, interfacial model. 

A more general formalism, described by COSP for any of its configurations, defines the mass 
spalling fraction, F s , as a function of the mass of oxide present before cooldown, W r ': 

F s = Q 0 W\ a (4) 

where Q 0 is a constant and the spall exponent a, usually fixed at 1.0, reflects the spalling 
dependence on scale thickness suggested by independent experiments [2]. Note that for the special 
case of a = 0, the spall fraction F s is a constant. 

The general sequence of the algorithm is that the scale is grown for a cycle duration of At, and a 
portion is allowed to spall according to eq. 4 above. The amount of retained scale is calculated and 
the growth upon the next cycle must be determined. For the latter, COSP uses the concept of 
effective time, t e ff, which is equivalent to that needed to grow the retained mass of scale in an 
isothermal exposure. Thus the total amount of new scale follows from the growth law selected from 
eq’s. 1-3, when teff+ At is now substituted for the time, t. 


The inputs, as called for in the publicly available software, COSP for Windows [4], thus 
require choosing a growth model and spalling configuration and specification of growth constants 
(k, m, and c), cycle duration (At), oxide phase or the stoichiometric constant (S c = mass oxide/mass 
oxygen in oxide), spall constants (Q 0 ) or area fraction (Fa), and spall exponent (a). For the Monte 
Carlo option, spall depth fraction (R spa u) and number of segments (n) are also required. Interfacial 
spallation is therefore handled by using an R spa ii=1.0. The primary output is the net sample weight 
(mass) change (AW/A) per cycle, but also includes the total amount of oxygen and metal reacted 
(£W 0X , EW m ), the current amount of retained scale (W r ) and the total amount and current fraction of 
scale spalled (XW S , F s ). No mechanistic assumptions are required (critical values of oxide 
thickness, strain energy, sulfur content, aluminum content, etc.). 



Figure 2. Trends with input parameters for COSP 
cyclic oxidation weight change curve [8]. Example 
shown for parabolic growth, uniform layer spalling 
for A1 2 0 3 : kp = 0.01 mg 2 /cm 4 hr, At = 1 hr. 

On = 0.002 cm 2 /mg, m = 2.0, a = 1.0. 



Figure 3. The effect of oxide type (S c ) on model cyclic 
oxidation curves [4]. Uniform layer spallation: 

(k p = 0.01mg 2 /cm 4 /hr, Q 0 = 0.002 cm 2 /mg, At=l hr). 


Cyclic Oxidation Responses. 

A typical model weight change vs 
number of cycles curve. Fig. 2, 
exhibits features similar to 
experimental curves: a maximum in 
the weight change, followed by a 
decrease to zero weight change and a 
final linear rate of loss. Regular 
trends in this and other outputs have 
been described in more detail 
previously in [1-5] and are suggested 
by the arrows. The maximum 
increases with the growth factor k p At, 
but decreases with metal content in 
the scale S c and spall constant Q 0 (or 
generally increased by all these 
parameters. The number of cycles to 
reach maximum (j ma x) and cross zero 
(j 0 ) are decreased by all these 
parameters for cases where the COSP 
spall exponent ooO, but remains 
invariant with k p At for the constant 
area fraction models [1,5], for the 
fixed probability models [6,7], and for 
COSP cases with a=0. It is found 
that the ratio of jo/jnm is exactly 3.0 
for most DICOSM cases and 
approaches 3.3 for the COSP uniform 
layer model (for AI2O3 and ot=l), thus 
in a sense determining the general 
form of the weight change curve. 

The actual response for various 
common oxides are shown in Fig. 3. 
Because the metal mass content 
increases with S c by definition, the 
effect on spalling is manifested here 
as higher rates of weight loss, even 
for the same growth rate, as in Fig.3. 


A more subtle effect is the decrease in the jo/jmax ratio, thus changing the ideal form of the curve 
as well. In fact, for the most general COSP models, it can be shown that this jo/jmax ratio, or shape 
index, is not really a constant at all. This effect is more explicitly presented in Fig. 4, showing 
significant decreases in the index with increasing stoichiometric mass factors, S c . Furthermore, it is 
seen that the growth exponent (linear, parabolic, or cubic with m=l, 2, or 3 in eqn.l) greatly 
increases the shape index. Finally, the spalling exponent, a, also has an effect. 

Thus low metal content scales, growing 
with sub-parabolic kinetics, and with a 
constant spall fraction can be expected to 
exhibit the highest shape index and a 
cyclic weight maximum greatly skewed 
towards the origin. To illustrate the shape 
factor for the most popular case, it can be 
seen that an S c of 2.124 for AI 2 O 3 , with 
parabolic growth and a spalling exponent 
of 1 , indeed produces a shape factor of 
3.27 as discussed above. While these 
basic variations may be disconcerting to a 
cyclic oxidation modeling effort, it is 
noted that the shape factor is nevertheless 
invariant with a number of other model 
parameters such as k. At, or Q 0 . 

Other outputs. The models also 
calculate various other elements germane 
to the cyclic process. The total amount of 
oxygen reacted, aluminum consumed, and 
oxide spalled are shown for the case of 
uniform spallation of an AI 2 O 3 scale in 
Fig. 5. They can be seen to reach linear 
steady state rates of increase at essentially 
the point where the weight change curve crosses zero. Also shown is the weight of the retained 
scale after spallation (W r ) and the fractional amount of spalling (F s ). These reach essentially 
constant values as part of the steady state regime. 

Comparison of various spalling conditions. From a systematic study of the various spalling 
formalisms over a range of input parameters, it is generally observed that the COSP uniform layer 
spalling model produces the least severe metal consumption rate for equivalent spall constants (Q 0 ) 
or with equivalent spall fractions (F s ). The DICOSM model can be most directly compared to 
COSP models where the value of the spall exponent a = 0, i.e., the spall fraction is a constant, eq. 4. 
Here, for example, it is found that the terminal steady state rate of metal consumption for the COSP 
uniform layer spalling is about 20% less than the interfacial bimodal spalling (Monte Carlo) case 
and about 15% less than the interfacial DICOSM case. For DICOSM and these special COSP 
cases, i.e., where a = 0 , it is also found that the number of cycles to reach maximum and zero 
weight change do not vary with the input growth parameters (k p . At). Basically, the curves are 
amplified vertically without changing the basic shape [2-5]. 



(S c -1), metal/oxygen weight ratio 


Figure 4. Variations in shape factor of COSP model 
weight change curves: decreases with stoichiometric 
factor, S c ; increases with growth exponent (m=l,2,3) 
for linear, parabolic, cubic kinetics; and decreases 
with spall exponents, a=l,2, 3. 



Figure 5. Various outputs of COSP cyclic oxidation model. Steady 
state behavior (i.e., limiting oxide thickness and linear rates of 
consumption and spalling) is achieved at ~ j 0 . (Same model parameters 
as in Fig. 2 [2-4]). 


DICOSM Equations and Characteristic Features 

In order to demonstrate more explicit functional relationships between cyclic oxidation responses 
and the input parameters, some attention will now be given to the DICOSM model formulae. In the 
model derivation, the outputs have been expressed in terms of series summations involving the 
square root of an integer [5]: 
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applying respectively to Case A, for cycle number j < n 0 , and Case B, for cycle number j > n 0 . 
Here n 0 is the number of effective area segments, only one of which spalls per cycle, and is 
equivalent to the inverse of the area spall fraction. Fa- Accurate, closed form solutions have been 
obtained through the use of an algebraic substitution for the summation series, ZjVi = Vz j + h j , 
the Good-Smialek Approximation (GSA) [5]. 

The resulting algebraic expressions may be solved for (AW/A) max , j max , jo, and the terminal slope 
of mass loss, T.S [5]: 


"AW" 

A /GSA, max 


Ay 


k At 
P 


(2(2S c -l))^ 


(7) 


J 


max 


fc-s c 

2 (2S C - 1) 


( 8 ) 


[ 9 ] 


j .. 3 (X~ S ° 

”' A = 2 (2S 0 -l) 

T.S. = -(S 0 -l\/F^at [10] 


These relations explicitly define how key features (descriptive or characteristic parameters) vary 
directly with the model input variables. For example, it is easy to discern that the number of cycles 
to reach zero weight change, j 0> is exactly 3.0 times the number to reach maximum, j max , and that 
each of these terms vary inversely with area spall fraction. Fa- The maximum in weight change 
varies with the square root of k p At and a more complex function of Fa and S c . Finally the steady 
state rate of weight loss per cycle varies with S c and the square root of the product of the spall 
fraction, growth rate, and cycle duration. Other attempts at expressions for characteristic 
parameters have been developed empirically and show some similarity to those above [1-7]. 

Universal (Normalized) DICOSM Cyclic Oxidation Behavior 

Weight change curves. The simplicity of these DICOSM-GSA equations suggest some special 
constructions to describe general cyclic oxidation behavior. Dimensionsless, universal terms for 
weight change, W u , and cycle number, J u , are obtained by normalizing by corresponding (AW/A) max 
and j^, yielding the cyclic oxidation curve for Case A as [5,8,9]: 

W„=X[3Jf-Jf] (ID 


This equation thus defines the 
universal cyclic oxidation weight 
change behavior for all values of the 
input parameters, up to the point 
where j > n 0 , as shown in Fig. 6. Here 
the normalized weight change equals 
zero at the normalized cycle number 
of 3, and the maximum occurs at J u = 
1. The branches for the steady state 
portion, Case B, begin at j = n 0 , which 
varies on the universal plot for 
differing oxides (shown here for S c = 
1 to 5). The steady state slopes must 
also vary (eqn. 10) and were verified 
empirically from running many 
combinations of input parameters. 


Material consumed. 

Now we consider the relations describing the amount of metal consumed, which is the pertinent 
quantity needed for failure criteria based on alloy depletion, mechanism changes, or breakaway 
oxidation [8,9]: 

zw - (s c - 1) f a jkm K i* + y 3 r } (12A> B) 

SW. G “ = (s o - 1 ) F a V^{in/ ! + X n ° ! ' ! } 



Figure 6. Universal DICOSM curves in which the weight 
change and cycle number are normalized by the coordinates 
at the maximum weight. A single curve is obtained for any 
kp, At, or F a for each S c specified [9]. 


These equations may now be used to represent universal metal consumption for all DICOSM 
model parameters. The cycle number is normalized by the number of segments, n 0 (=1/Fa), and the 
normalized weight change, W met ,u, is defined as eq. 13. The normalized metal consumption 
relations for Case A and B (eq. 14A, B) are produced, giving the single, continuous plot in Fig. 7: 
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In essence, this construction defines the amount of metal consumed for any combination of S c , k p . 
At, and F A , with the transition to steady state behavior defined by the coordinates (1,1). 
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Normalized cycle number, J u = j / n 0 = jF A 


Figure 7. Universal, normalized 
curve for metal consumption. 
Universal transition point at the 
coordinates (1,1) and log-log 
slopes of 1/2 before and 1.0 
after the transition [8,9]. 


Failure mapping. Finally, a universal construction has been applied to the concept of cyclic 
oxidation failure maps, modified from those first presented by Poquillion and Monceau [6,7]. Now 
a (critical) prescribed lifetime and prescribed amount of material consumed is presented on a spall 
constant, growth rate plot, or F A -k p diagram in the case of DICOSM. Again, by suitable 
normalization, a universal construction can be developed which defines the universal failure locus, 
once the specific failure criteria (j*, £W m *) have been selected [9]. Four regimes or pseudo- 
quadrants of behavior are distinguished, corresponding to survival (I), or failure controlled by 
various combinations of spallation and scale growth (n,n,IV). By plotting published experimental 
data on such a map, for a failure criterion of 10 mg/cm 2 metal consumption in 1000 hr, it is seen 
that using the cyclic 1200°C oxidation data, MoSi 2 [10] and Zr-doped NiAl survive, whereas 
undoped NiAl fails primarily by spallation [11]. Similarly, the effect of sulfur on increased scale 
spallation had been plotted as a trend toward higher spall probabilities on a p-k p map [7]. Finally, 
for a less demanding 100 hr life criterion, pure nickel is seen to fail by excessive growth at 1200°C 
[11] and pure copper fails at 500°C by contributions from both growth and spallation [12]. 


\ 



normalized k p , 
{j*k p At(S c -1) 2 /(3/4SW* m ) 2 } 


Figure 8. Universal F A -kp life map 
for DICOSM model, showing four 
major sectors and experimental 
results for diverse materials [9]. 

Survival appears as Region I, 
failure by spallation (II), growth 
(IV), or both (HI). 

Oxidative failure criteria set as 10 
mg/cm 2 metal consumed in 1000 hr. 
for MoSi 2 and NiAl(±Zr) at 1200°C, 
in 100 hr for Ni at 1200°C;and in 100 
hr for Cu at 500°C. 


Concluding Remarks 

This paper has reviewed the basic elements, assumptions, and outputs of cyclic oxidation models. 
Inherent in the COSP program is the flexibility to specify various scale growth kinetic laws, spall 
morphology, and thickness dependence. More severe material loss is seen to occur for interfacial 
spallation and when the spall fraction increases with scale thickness (a>0). A special deterministic 
model (DICOSM) enabled simple mathematical approximations and normalization by performance 
terms. This allowed the construction of universal cyclic oxidation weight change, total metal 
consumption, and failure mapping plots for a wide variety of input parameters, indicating their 
generalized forms. 
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